f(R) Gravity and Chameleon Theories 
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We analyse f(R) modifications of Einstein's gravity as dark energy models in the light of their 
connection with chameleon theories. Formulated as scalar-tensor theories, the f(R) theories imply 
the existence of a strong coupling of the scalar field to matter. This would violate all experimental 
gravitational tests on deviations from Newton's law. Fortunately, the existence of a matter depen- 
dent mass and a thin shell effect allows one to alleviate these constraints. The thin shell condition 
also implies strong restrictions on the cosmological dynamics of the f(R) theories. As a consequence, 
we find that the equation of state of dark energy is constrained to be extremely close to —1 in the 
recent past. We also examine the potential effects of f(R) theories in the context of the Eot-wash 
experiments. We show that the requirement of a thin shell for the test bodies is not enough to 
guarantee a null result on deviations from Newton's law. As long as dark energy accounts for a 
sizeable fraction of the total energy density of the Universe, the constraints which we deduce also 
forbid any measurable deviation of the dark energy equation of state from -1. All in all, we find 
that both cosmological and laboratory tests imply that f(R) models are almost coincident with a 
ACDM model at the background level. 
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I. INTRODUCTION 



The acceleration of the Universe expansion was dis- 
covered ten years ago and is still a deep mystery (see e.g. 
[l[ for recent results on observations of dark energy and 
e.g. [2j, [3J for theoretical overviews). Two types of ap- 
proaches have been considered. One can either introduce 
a new kind of matter whose role is to trigger accelera- 
tion or modify the behaviour of gravity at cosmological 
distances. In the first approach, dark energy is a new 
energy form, with all its well-known puzzles such as the 
cosmological constant problem, the coincidence problem 
and the value of the equation of state. This approach is 
subject of intense experimental investigation and any de- 
viation from -1 would be a smoking gun for new physics 
beyond the standard models of particle physics and cos- 
mology. On the other hand, in the second approach, 
various attempts to modify gravity have been presented 
(see e.g. the literature is vast, see for a re- 

cent overview and further references). Up to now, they 
are all plagued with various theoretical problems such as 
the existence of ghosts or instabilities. In this paper, we 
will consider a modification of Einstein's gravity, the so- 
called f(R) theories, which do not seem to introduce any 
new type of matter and can lead to late time accelera- 
tion. In fact, these theories can be reformulated in terms 
of scalar-tensor theories with a fixed coupling of the extra 
scalar degree of freedom to matter. As theories of dark 
energy, they suffer from the usual problems and are also 
potentially ruled out by gravitational tests of Newton's 



law. 



The only way-out for these models is to behave as 



chameleon theories 116 
pendent mass [TtI E 



i.e. develop an environment de- 
ij, [20| . When the density of the 
ambient matter in which the scalar field/chameleon prop- 
agates is large enough, its mass becomes large and the 
smallness of the generated fifth force range is below the 
detectability level of gravitational experiments. On the 
other hand, for planetary orbits or any other situations 
in which gravity is at play in a sparse environment, one 
must impose the existence of the so-called thin shell ef- 
fect. In this case, the fifth force is attenuated as the 
chameleon is trapped inside very massive bodies (the 
Sun for instance). It has been argued that the existence 
of thin shells is usually enough to salvage f(R) models 
Ojj]. We show that thin shells do not always guar- 
antee null results in experimental tests of Newton's law. 
We exemplify this fact using the Eot-wash setting and 
obtain strong constraints on the models which translate 
into stringent bounds on the present dark energy equa- 
tion of state, preventing any detection of a deviation from 
-1 in the foreseeable future |l+io| < 1CU 4 , where w is the 
equation of state of dark energy in the recent past. This 
corroborates a similar bound obtained from the existence 
of thin shell for objects embedded in a super-cluster. It 
should be noted that this result holds at the background 
level. For higher red-shifts where the effective dark en- 
ergy density fraction, f2de, may become small (or even 
vanish), larger deviations can be present as exemplified 
in the models in [2l|, [52] where the equation of state can 
deviate from -1 for red-shifts of order z = 2 — 3. In all 
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these models however |1 + u>|r2dc <C 1, and so even if 
w deviates significantly from — 1, deviations of the ho- 
mogeneous cosmology from ACDM are still very small. 
Detectable deviations from ACDM are envisageable at 
the perturbative level as the growth factor is anomalous 
at small scales (see e.g. |23| for a discussion of this point 
for the original chameleon model). Some consequences 
of this fact on the matter power spectrum and the CMB 
spectrum of f(R) models have been presented in Ref. 

The paper is organized as follows: In the subsequent 
section, we review f(R) models and chameleon theories. 
In section III we derive the cosmological thin shell bound 
on the equation of state. In section IV, we consider tests 
of the inverse square law. Finally, we apply these consid- 
erations to particular models in section V. The appen- 
dices contain some technical details. 



II. f(R) GRAVITIES AND CHAMELEON 
THEORIES 



A. f(R) theories 

An f(R) theory is a modified gravity theory in which 
the usual Einstein- Hilbcrt Lagrangian for General Rela- 
tivity, i.e. R, is replaced by some arbitrary function of 
the scalar curvature i.e. f(R). The action for an f(R) 
gravity theory therefore takes the following form: 



= d 



9^ i f(R) + Sm a tter\9 f ,u,^i], (1) 



where the \tj represent the matter fields. 

In this work we are concerned only with metric f(R) 
theories, in which only the metric g^ v is the independent 
variable in the gravitational sector. The quantity Y p v is 
taken to be the Levi-Civita connection associated with 
the metric g^. In these metric f(R) theories the field 
equations are: 

R^f(R) -\f{R) 9liV = KT™^ (2) 



-V i y v f{R)-g ttv Uf{R). 



where k = 1/M, 



B. Transformation to a Scalar- Tensor theory 

Eq. © gives a set of equations which are second or- 
der in derivatives of R, which is itself second order in 
derivatives of g^, making the field equations fourth or- 
der in g^„. Finding solutions to fourth order equations 
can be mathematically and physically troublesome, but 
fortunately metric f(R) theories can be recast as a scalar 
tensor theory with only second order equations via a well 
known conformal transformation. We define <f> by 



exp 



Mpi 



) = f'(R), 



where (3 = \/l/6. We also define the Einstein frame 
metric g M „ by a conformal transformation 



and let R be the scalar curvature of g^ u . When rewritten 
in terms of g^ and <f>, Eq. ((!]) becomes: 



Sst = 




M 2 _ 1 



matter [e M ^ g^^i], 

where the potential V((j>) is given by: 



V(<f>) = 



M^(Rf'(R)-f(R)) 
2f'(Rf 



(3) 



(4) 



When the action is written in the form of Eq. ©, we 
say that we are working in the Einstein frame. The field 
equations then become: 



= Rfj, v - -Rg^v = kV^V^ 



□0 = v'(<p) 



/jimatter 



JL 

Mpi 



(5) 



(6) 



In the above and subsequent expressions, the covariant 
derivatives, V M , obey V^g^^o and all indices are raised 
and lowered with g^ LV unless stated otherwise. We note 
that in the Einstein frame T'™ atter is not conserved but 
instead: 



VT-imattcr fi 



JL 

M P f 



rpTi 



(7) 



This implies that matter will generally feel a new or 'fifth' 
force due to gradients in <f). We note from Eq. (0) that, 
when written as a scalar tensor theory, gravity in an f(R) 
theory is essentially General Relativity, and all the mod- 
ifications are essentially due to the effective 'fifth force' 
and to the energy density of 4>. Much of our intuition 
for how gravity works is based on how it works in Gen- 
eral Relativity. When an f(R) theory is written as a 
scalar tensor theory we can readily make use of this intu- 
ition in solving the field equations. This may not be the 
case, however, in the original frame in which the equa- 
tions were fourth order and so in those circumstances one 
would have to be more careful. Note that all physical ob- 
servables must be independent of the choice of frame, i.e. 
the choice of metric g^ or g^ v . 



C. Chameleon Theories 

Since f(R) theories are equivalent to scalar tensor the- 
ories one can generally directly apply the plethora of con- 
straints on scalar tensor models. In particular, since the 
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extra degree of freedom, <f>, couples to matter with grav- 
itational strength, tests of the inverse square law such as 
the Eot-Wash experiment require that <j) have a mass, 
= ^jV^cf,, greater than 1 meV. Cosmologically 4> 
would then have been fixed at its minimum since very 
early times, and physics over astrophysical scales would 
be indistinguishable from predicted by unmodified Gen- 
eral Relativity with a cosmological constant. Both the 
coincidence problem and the problem of the small size of 
the cosmological constant would not be alleviated in this 
scenario. However, this is not the whole story. Labora- 
tory constraints on scalar tensor theories can be greatly 
relaxed if m,p = \JV^ develops a strong dependence on 
the ambient density of matter. Theories in which such 
a dependence is realized are said to have a chameleon 
mechanism and to be chameleon theories. In such the- 
ories, <j) can be heavy enough in the environment of the 
laboratory tests so as to evade them, whilst remaining rel- 
atively light on cosmological scales. It must be stressed 
that even with a chameleon mechanism, it is still very 
difficult, if not impossible, to construct such a theory 
where the late time cosmology would be observational 
distinguishable from the usual ACDM model. To the 
best of our knowledge all such theories which are also 
experimentally viable require a fairly high degree of fine 
tuning to ensure that the effective cosmological constant 
is small enough. 

A chameleon theory is essentially just a scalar-tensor 
theory in which the potential has certain properties. As 
such Eqs. ([3] -[7]) also define a chameleon theory for cer- 
tain classes of V(</>). In these circumstances the f(R) 
theory would be equivalent to a chameleon theory. In 
a general chameleon theory, /3, which parametrizes the 
strength of the coupling of <f> to matter, could take any 
value and potentially even be different for different mat- 
ter species. If a chameleon theory is equivalent to a f(R) 
theory, however, f3 is fixed to be y/l/6 and is the same 
for all types of matter. If a f(R) theory is not equiva- 
lent to a chameleon theory it would be generally ruled 
out by laboratory tests of gravity and / or result in no 
detectable deviations from General Relativity over astro- 
physical scales. 

For an f(R) theory to have a chameleon mechanism 
one must require that, in at least some region of </>-space: 

V'{<t>) < 0, V"(cf>) > 0, V'"{<p) < 0. 



It follows from the definition of d> that: 



dR 



Mpi f" 



2/3 /' 

and therefore the derivatives follow 
0M Pl 



V(4>) = 
v"W = 

V"'W) = 



r 



[Rf - 2/] 



R 
J' 



3Mpi 



1 

3 



f f" 

-J7T 



R 

7 



8/ 
f' 2 



(8) 

(9) 
(10) 

(11) 



In general, this gives strong constraints on the form of 
f(R)- In the following we will study examples where 
these conditions are met. 

When these conditions are satisfied, the mass of (j> in 
a suitable large region with density p will increase with 
p. In order to evade constraints coming from local tests 
of gravity, it is not, however, enough that a theory pos- 
sess a chameleon mechanism; the mechanism must, in 
addition be strong enough for chameleonic behaviour to 
occur for the test masses used in the laboratory gravity 
experiments. 



D. Thin-Shells 

1. Chameleon Theories 

Chameleon theories do not behave like linear theories 
of massive scalar fields. In situations where massive bod- 
ies are involved, the chameleon field is trapped inside 
such bodies and its influence on other bodies is only due 
to a thin shell at the outer edge of a massive body [HI] . 
As a result, the field outside the massive body for dis- 
tances less than the range of the chameleon force in the 
outer vacuum is effectively damped leading to a shielded 
fifth force which becomes undetectable. The criterion for 
a thin shell is 



Ac/) 
mpi 



< 2/3$ 



(12) 



where A<f) — (f)^ — 4>o 1S the field difference from far inside 
the body to very far away. We define the body and the 
region outside it to have densities po and p^ respectively. 
It involves Newton's potential $n &t the surface of the 
body. In general the field values at infinity, </>oo , and deep 
inside, 4>o, are related to and po by 



dsV 



-/ML. 
m P i 



(13) 



In most current situations involving runaway potentials, 
when po ^> poo, this implies that tfioo > (po- Hence, 
A</> = 4>oo implying that cosmological information can be 
inferred from local tests. Moreover, in a cosmological set- 
ting, the chameleon sits at the minimum (|13[) during the 
matter era. As a result, the variation of the equation of 
state in the recent past is severely constrained. Another 
important consequence of the chameleon effect is the ex- 
istence of an anomalous growth of the density contrast 
for scales lower than the inverse mass of the chameleon, 

•4. Tl V V, ~ -l+- v /l+24(l+2 / 3 2 ) ir^l T 

i.e. it grows like a where v w — = 123JJ- In. 

the f(R) setting, some of the consequences of this anoma- 
lous growth on the CMB and the matter power spectrum 
have been analysed using the convenient variable 



B 



[rr dR 
fn dH J 



H 



(14) 



whose square root represents the compton wave-length, 
i.e. the inverse mass of the chameleon, in horizon units 
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i? _1 [24|, [H, HI]. Effects on structure formation could 
be seen for values as low as B = 10~ 4 in future galaxy 
surveys [IB} . In the following we will find an explicit ex- 
ample of logarithmic f(R) model which could lead to 
effects on scales as large as 100 h~ 1 Mpc. All these facts 
will be crucial in the following. 



2. Thin shells in the language of f(R) theories 

It is useful to write the function f(R) in the form 
f(R) = R + h(R), where h measures the deviation to 
Einstein gravity. To leading order, as a consequence of 
Ref. [l6|, the thin shell condition can be formulated as 



\Ah'(R)\ < -<J> W . 



(15) 



As Newton's potential is small on cosmological scales, 
with an upper bound around 10~ 4 , this implies that h' 
must have very small variations. The thin shell condition 
is a constraint on local experiments at the present time. 
It has nothing to say, a priori, about the evolution of 
the universe since matter equality for instance. Another 
useful combination (which is not to be confused with the 
chameleon mass m^) has been used 



Rti'(R) 
1 + h'(R) 



(16) 



It has been shown that the existence of a matter era fol- 
lowed by an accelerated period requires m < 0.1. For 
models where m is (nearly) a power law, the thin shell 
constrain implies that m is much smaller for reasonable 
powers. In the following, we will obtain a bound on the 
equation of state at present time which implies that de- 
partures from ACDM are tiny. 



III. THIN-SHELL CONSTRAINTS ON 
COSMOLOGY 

In subsequent sections we will assume that test bod- 
ies used in laboratory based gravity experiments have 
thin-shells. In the absence of any thin-shell, the inverse 
square law tests, such as the Eot-Wash experiment [32[ 
(as well as other tests of gravity over longer ranges), rule 
out theories with f3 = 1/^/6 as it is in f(R) theories. The 
thin-shell requirement must therefore be satisfied by any 
physically viable f(R) theory. Although it is not often 
appreciated, the thin-shell condition for laboratory test 
masses actually places extremely tight constraints on the 
recent cosmological evolution of <fi. In this section we con- 
sider those constraints in the context of a general f(R) 
theory. 

In any single field scalar tensor theory there is a choice 
of frame. In the Jordan frame, the laws of physics in 
a local inertial frame are the same everywhere, however 



Newton's constant, Gjv, is different at different points in 
space and time. In the Einstein frame, Gn is chosen to 
be fixed but, as a result, local particle physics is position 
dependent. The process of converting astronomical ob- 
servations to cosmological parameters generally involves 
making assumptions about how today's laws of particle 
physics are related to those in the past. This said, if the 
relative changes in Gn (in the Jordan frame) are small 
i.e. <C 1, the differences between cosmological parame- 
ters in the two frames are only very slight. For instance, 
to calculate a redshift, one must compare the observed 
wavelength, A b s of a particular absorption or emission 
band to the wavelength that band would have had at 
emission, A c . Since one cannot go to the astronomical 
object in question and directly observe the wavelength 
at emission, it is generally assumed that particles physics 
in the past obeyed the same laws as it does today and 
so replace A c with the wavelength of the band as it is 
measured in a laboratory today, A to day When one is 
dealing with scalar-tensor theories, the assumption that 
\ = \oday is equivalent to a choice of frame, specifically 
the Jordan frame. 

To make comparison with observations straight- 
forward, one should therefore quote cosmological param- 
eters for the Jordan frame. This said, it is often more 
straightforward to perform calculations in the Einstein 
frame and then merely express the results in terms of 
Jordan frame quantities. 

Cosmologically, in the Jordan frame we have: 



ds 2 



fl2 ( 7 ?) [ — d?? 2 + TijdxMa;- 7 '] 



(17) 



and <f> obeys: 



1 r 
—^Q.rin ~ 2— hr^.ri = 77 T ma ttor 

- a 6 3 



a 



+ 2<F 3 U$] (18) 



where $ = e ~ 20 ^ Mpl = f'(R). At late times, when it 
is appropriate to ignore the contribution of radiation to 
the total energy density of the Universe, we have 

3^ = K^tter _ 3o^ _ 3fc 

a 4 $ a 3 <E> a 2 

The Einstein equations also give: 

_ < = kW(0) _ |m + ™ (20) 

a 3 a 4 <I>a 2 a 2 a 3 3> 

We assume that measurements are interpreted in terms 
of General Relativity, where the energy density of the 
Universe is assumed to be due to non-interacting, dark 
energy and normal matter. Thus we write 



fj2 _ a ,n _ K . K eff / 9 -n 

H ~& 3$/ mattcr + 3$/ dc " ^ (21) 



The above equation partly defines p*^, fi^f and Oj*; to- 
day $ = $o- Now in the Jordan frame: p ma ttor oc a~ 3 . 
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If the effective dark energy equation of state parameter, 
iu d §, were constant it would obey: p 1 ^ cx a~ 3 ( 1+w *<<\ 
More generally however the effective dark energy equa- 
tion of state is then given by: 



Taking the 77-derivative of Eq. (|2T|) we get: 

•,2 



(22) 



3$o 



Pmatter,77 



and so using the Eq. 



and Pmatter oc a 3 we have 



— T~ 2 ~ "^-Pmatter - ^-(1 + w dc)Pdc- 

gt a 4 a z <3>o "Po 

Finally by adding 3H 2 to both sides and using Eq. (|2~lj) 
we have: 

a % , & _ K „„eff nff 

So by rearranging the Friedmann equations we have 
found that 



eff „ „eff , 







a 3 


a 4 






$a 2 


$a 3 



(23) 



- K$y(</>) 



By comparing Eqs. (|19p and (f2"Tj) we see that: 



-« eff 



$ ac $ 
Therefore, 

(l+ W f)^f/$o 



/<f>0 iU K $y($)-3- fl ^- 



"/'matter 



a 3 $ 



$a 2 



2$, r^.r? 
$a 3 



(24) 



K ( $0 

+ Pmatter - r 1 

$0 V $ 



Thus, using 3f7^i? 2 = Kp e £/$ , we have: 



(l + wf) nf c 



'/<; 



2$. 



3$a 2 ff 2 3$aff 
*° -1 ' ()-»■ 



(25) 



(l + ) 17^ parametrizes the magnitude of deviations 
from ACDM. If = -(M P i/2/3) ln$ has changed by A</> 
in the last Hubble time, Eq. (|2l)|) implies that, in the 
recent past and in the Jordan frame, to within an order 
of magnitude: 



For later use we rewrite Eq. (|25p in terms of p = In a 



2*„ 



cff\ QCff 

w dc ) J 'de 



elf 



3(2 



(f - 1) [g 

2 + % 



(27) 



In both the Einstein and Jordan frames, e m pi — 1 
gives the relative change in the ratio of any particle mass, 
and the Planck mass, Mp\ between the times when 
4>i and when 4> = 4> . In the Einstein frame Mpj 
is constant but m p varies whereas in the Jordan frame 
the converse holds; the ratio of the two masses, being 
a dimensionless quantity, is the same in either frame. 
WMAP constrains any such variation in m p /Mp\ be- 
tween now and the epoch of recombination to be < 5% 
at 2cr (< 23% at 4a) (35[. It follows that since recombi- 
nation 



|e M pi - 1 < 0.05. 



(28) 



Light clement abundances provide similar constraints 
on any variation in Newton's constant between the 
present day and the time of nucleosynthesis [36| . 

Thin-shell constraints, however, provide an even 
tighter bound on the allowed change in (f>. To consider 
these constraints we work in the Einstein frame, however 
A(f> is the same in either frame. 

We assume, as is the case for the real Universe, that the 
scales of the inhomogeneous regions are small compared 
to the horizon scale, and that the Universe is approx- 
imately homogeneous when coarse-grained over scales 
larger than some L^ om <C if -1 . Thus over scales larger 
than I/hom, 4> ~ 4>b(t) and since Lhom -C H~ l , we can 
work entirely over sub-horizon scales, which simplifies the 
analysis greatly. We also assume that the curvature of 
spacetime is weak over scales smaller than Lh om - This is 
equivalent to assuming that the Newtonian potential, U, 
is small as are the peculiar velocities, v l , of any matter 
particles, i.e they are non-relativistic. 

Exploiting both the assumption that i/Lhom <C 1 and 
that gravity is weak inside the the inhomogeneous regions 
i.e. U <C 1 and v l v l < 1, we write <j) = (pb(t) + 5(f) and 
have to leading order in the small quantities and over 
sub-horizon scales: 

V 2 S<j) = V'(4>) + ^f- + 4>b + 3H<j> b . 



Now 



and so 



Mpi 



-4> b -3Hj> b = V, 4> {0 b ) + 



Mi 



pi 



V Mpi J 



(26) 



V 2 ^ = m 2 <50 + + A(<f>, <t> b ). 



Mpi 



6 



where 



Thus 



A{4>,<t>b)= [V'^b 



- VI - m 2 b t 



i r „-m\x-x'\ 



Air J \x — x' 

-A{ct>{x',t),<t> b {t))] 



f3Sp(x',t) 



M, 



pi 



(29) 



It is straightforward to show that the condition V'" < 0, 
which must hold for any chameleon theory, implies that 
A{(f), (f> b ) < for all <f> and 4> b . Thus 



> 



1 { aS ^ 

-— / dV — 

An I \x 



l^'l f36p(x',t) 
~\ Mpi ' 



Now if we require that a test mass at r = with central 
density p c > p b has a thin-shell, we must impose that at 
r = 0, (f> w cj) c , where 



V, 







Mi 



pi 



Thus must be able to change by at least <j) c — 4> b — 
—A(j) bc < i.e. we have the following necessary condition 
for thin-shell 



f3Acj> bc 



M, 



< 



pi 



Mi 



pi 



p — m\x — x | 

dV- — GSp[x',i). 



x 



x' 



(30) 



The right hand side of this equation is 0(U/3) or smaller, 
and the largest values of the peculiar Newtonian potential 
for realistic models of our Universe are roughly < 10 ~ 4 , 
and are generally around 10~ 6 — 10~ 5 for large clusters 
and superclusters [26j |. Thus we have the following con- 
servative constraint on the cosmological value of the field 
today: 



be 



M, 



< 10" 



pi 



(31) 



We have defined f(R) = R + h(R). The thin- 
shell constraint certainly ensures that cosmologically to- 
day \/3(j)/Mpi\ <C 1 and since we have 1 + h'(R) = 
exp(— 2/30/Mpi) by definition we are therefore justified 
in assuming that we have |ft'(i2)| <C 1- Then assuming 



that \h' 
by: 



and 



< 1 wc find that the potential, V{(f)), is given 



KV(<f>)n-(Rh'(R)-h(R)), 



1 



-V, 



R(l - 2h\R)) + 2h(R) - Rti(R). 



0Mpi 

To leading order then in \h'(R)\ we have: 



The chameleon mass squared, rr& = is then given, 
to leading order, by: 



= V, 



3h"(R)' 



(32) 



Provided m 2 ^/H 2 <§; 1, then the chameleon field will re- 
main close to the minimum of its effective potential [29l | 
cosmological, i.e. = — /3p m atter/Mpi and the energy 
density of the chameleon field will be dominated by its 
potential. Assuming that this is the case we would have: 

R « AnV{4>) + 

and defining ft m = K,p m&ttcv /iH 2 and Sl dc ~ kV(4>)/3H 2 , 
we have: 

R^3(AQ dc + n m )H 2 , 
and so m 2 ,/H 2 ^> 1 becomes: 



■fin 



Rh"(R) 



> 1. 



Therefore, in many theories, an observationally viable 
evolution of 4> requires that it has sat close to the effective 
minimum of its potential since recombination [29l | i.e.: 



/3pmatter(i) 



M, 



PI 



Since the background density of matter decreases with 
time, V^cj, > implies that 4> increases with time. Thus 
for test mass with density p c ~ 0(1) gem -3 , we have in 
the recent past, i.e. out to z « 1: 

(j) c < <p b (t) < <p b (t ) 

where t = to is the current time. In this case Eq. (|3"Tj) 
gives the following conservative constraint: 



(&(*o) ~ &(*)) < 10" 



JL 

M P , 



and so, from Eq. (|26j) we obtain that 
\l + w efl \Qf e < 10" 



(33) 



In the recent past where is not negligible, this leads 
to a stringent constraint on the deviation of the equa- 
tion of state from -1. It should be noted that although 
1 1 + w c ff | £1°^ is constructed simply out of the scale fac- 
tor, a, and its derivatives, neither w c s(z) nor f2^(z) are 
uniquely defined as functions of redshift in models such as 
these where the scalar field interacts with normal matter. 
As a result, it is possible to define fij? so that it vanishes 
and even becomes negative in the past. If such a defi- 
nition is made, then one would (unless the |1 4- w e g \ 
also happens to vanish) predict that w e g diverges, and 
hence deviates significantly from —1. A behaviour such 
as this was noted in Refs. [2l], [22J. As a result, an ap- 
parent effective deviation from ACDM can be deduced. 
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However, because of the freedom to redefine 0,^ and 
hence w e s, one should not rush to assign any physical 
meaning to the divergence of w e s, and deduce that it 
represents a significant deviation from ACDM, since one 
could always remove this divergence by choosing to de- 
fine fijjp it such a way that it is positive definite. In 
all cases, the bound (|33[) gives an intrinsic measure of 
the deviation of the background cosmology from ACDM, 
and it all cases it is small. Therefore, the predicted late 
time cosmology is observationally very close to ACDM. 
Additionally, the prospects for being able to detect such 
small deviations for ACDM at the background level in 
the near future are poor. Of course, as we have already 
mentioned, detectably at the perturbative level might be 
within reach. 

This said, the thin-shell constraints do not themselves 
rule out larger deviations from ACDM. It may be that 
(3cf) /Mp\ has undergone relatively large changes in the 
past i.e. much larger than 0(1O -4 ), but that we now just 
happen to live at a point in time when f3A<fib c /Mp\ < 
1CP 4 . This would, however, be a fairly remarkable co- 
incidence and would inevitably require a great deal of 
fine-tuning of the theory and the initial conditions. To 
avoid this new coincidence problem, we would have to re- 
quire that the cosmological changes in /30/Mpi have been 
smaller than 0(1CP 4 ) in the recent past which would 
in turn, as we illustrated above, constrain any devia- 
tions from ACDM to be unobservably small. We note, 
however, that deviations can be expected on very small 
scales, as in the original chameleon model (23[. 

In this section we have sketched how the thin-shell re- 
quirement for laboratory test masses place a very strong 
constraint on the recent cosmological evolution of <f>, and 
generally constrains any deviations from ACDM in the 
predicted cosmology to be small. This is not, however, a 
'water tight' constraint as it may be possible to circum- 
vent it by requiring a seemingly improbable cosmological 
evolution wherein such bodies would only have developed 
thin-shells in the recent (in the cosmological sense) past. 
The laboratory constraints which we will derive in what 
follows cannot be avoided in this way. 



scales r ~ l/m^. If, as in chameleon theories, is 
not a constant but instead undergoes 0(1) or greater 
variations, the behaviour of the force is more complicated 
but generally not of inverse square law form. 

It is often assumed that what is needed for an f(R) the- 
ory to avoid the constraints of inverse square law tests, is 
that the test bodies develop thin-shells. Generally, how- 
ever, this is not the case. The presence of a thin-shell 
causes the chameleonic force due to a body to drop off 
much faster than 1 /r 2 near the surface of the body. Far 
from the body, the force has a Yukawa form, although 
as a result of the fast drop-off near the surface, it is 
much smaller than one would normally expect it to be. 
If two thin-shelled bodies are sufficiently close however 
then they would be inside the region where the faster 
drop-off is occurring. In these cases the detectable vio- 
lation of the inverse square law can be much larger than 
one might expect. 

A number of different experiments have searched for 
violations of the inverse square law. For gravitational 
strength forces, i.e. a ~ 0(1), the best constraints are 
currently provided by the Eot-Wash experiment [32| . 

The Eot-Wash experiment [HJ consists of two plates: 
the attractor and the detector. The detector is 0.997 mm 
thick and made out of molybdenum. The detector has 42 
4.767 mm diameter holes bored into it in a pattern with 
21-fold azimuthal symmetry. The attractor is similar and 
consists of a 0.997 mm thick molybdenum plate with 42 
3.178 mm diameter, arranged in a pattern with 21-fold 
azimuthal symmetry, mounted on a thicker tantalum disc 
with 42 holes, each with diameter 6.352 mm. The holes in 
the lower tantalum ring are displaced so that the torque 
on the detector due to the attractor from forces, such as 
Newtonian gravity, that have a 1/r 2 behaviour vanishes. 
The detection of a non-zero torque would therefore indi- 
cate the presence of either a correction to gravity with 
a behaviour different from 1 /r 2 or the presence of a new 
force that also did not behave as 1/r 2 . 



IV. INVERSE SQUARE LAW CONSTRAINTS 

In the weak field limit, the gravitational force due to a 
small body drops off as 1/r 2 , where r is the distance to 
the body's centre of mass. If there is an additional scalar 
degree of freedom to gravity with constant mass m^, the 
force instead drops off as: 

(1 + a(l +TO r)e- m * r )) 

where a parametrizes the strength with which the scalar 
degree of freedom couples to matter. In f(R) theories 
a = 2f3 2 = 1/3. When m^r « 1 or m^r > 1, the 
force still drops off, approximately, as 1/r 2 , however there 
would be a noticeable deviation from this behaviour over 



A. Chameleonic Force & Torque 

We now calculate the force, due to a chameleonic scalar 
field, <j>, on one plate due to the other lying parallel to it. 
From this we calculate the chameleonic contribution to 
the torque. 

In a background, where <f> = (fib far from the plates, 
the chameleonic force per unit area between two parallel 
plates, of the same or similar compositions, both with 
thin-shells and with a distance of separation d between 
their two facing surfaces was found, under certain condi- 
tions, in Ref. [28|,[3l|]. In Appendix lAl we generalise those 
formulae. We find that the chameleonic force between 
two parallel circular plates, with radius r p and thickness 
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D <C r p , and separation d <C r p is given by: 



5= 



(V/ - V') 2 



(34) 



2m\C{m h r p ) 



where 0o (^) is defined to the values of <j) midway between 
the two plates, and formulae for it are provided in Ref. 
13111 . We have also defined: 



E = 
C = 

D = 



1 + 2C(m b r p )D((j) c , <f> b ) - Vl + 4CZ?, (35) 



(e 



m b r p 



1/2) 



m b r p 



(e 

< [v(<M - ~ W<0(& - &)] 
(^OW-Wc)) 2 



(36) 



(37) 



The last term in Eq. (|34|) represents the only difference 
between the generalized force formula and the one pre- 
sented in Refs. (HH) & ([ll|), and we note that the 
extra term is independent of the separation d. When 
C(rribr p )D((j) c , <j>b) <C 1, the last term in Eq. (j3~4"|) is neg- 
ligible. We note that C(mbr p ) -C 1 when rribr p ^> 1, and 
so whenever m,br p 3> 1, the last term is always negligible. 

The details of how F$(d)/A drops off with d will de- 
pend on the form of V ((/>). For many choices of V(<j>), 
e.g. V{4>) oc <j) n for n < —2 or n > 2, one finds that 
F<j,(d)/A drops off faster than 1/d, for all d > c?2 is small 
compared to both r p and the radius, r/j, of holes in the 
plates. Indeed, this will certainly be the case, no matter 
what form V(<fi) takes, if m^r/j ^> 1 where m b = m^ipb)- 
Provided this is the case, we can define the potential en- 
ergy, V$(d) due to the chameleonic force for two plates 
with separation d -C Rh thus: 



V^(d) w A 



A 



(38) 



The faster 1 jd drop off has been used to set the upper 
limit of the above integral to oo. 

In the Eot-Wash experiment the plates have a number 
of holes in them. This means that as one plate is rotated, 
by an angle 9 say, relative to another, the surface area, 
A((f>), of one plate that faces the other changes. Note 
that F^/A does not depend on A. The torque due to the 
chameleonic force is given by the rate of change of the 
potential V$(d) with 9: 



We therefore have: 



r° 

T^d) = a T 

Jd 



A 



A 



ds. 



ds, 



(39) 



(40) 



where ar = dA/d9 is a constant that depends only on the 
details of the experimental set-up rather than the theory 



being tested. For the 2006 Eot-Wash experiment [3S 
find 



we 



3.0 x 10" 3 m 2 . 



If F^ I A drops off too slowly over scales of the order of 
rh then a more complicated analysis must be performed, 
and knowing the force between two infinite parallel plates 
is no longer enough to find a good approximation to the 
torque. Instead a full numerical analysis would have to be 
undertaken to get accurate results. This said, for d > r^, 
we do not expect F^(d) to depend strongly on 9 because 
the effect of the holes will be largely smeared out over 
separation distances much larger rh. On scales <^ r^, we 
found that F oc A{<j>). Since T = dV< t> /d9, dV^/dd = 
F^, and we expect F^ to be largely 9 independent for 
d ^> rh and oc A{9) on smaller scales, we expect, to 
within an order of magnitude, that: 



T$(d) w a T 



A 



ds. 



in these cases, where once again ay = dA/d9. By pick- 
ing rh as an upper bound for the integral we are probably 
under estimating the torque as we are dropping the con- 
tributions from larger separations. 



B. The Effect of an Electrostatic Shield 

Up to now we have not considered the role played by 
the electrostatic shield. Because the shield is so thin 
(d s = 10 /im) compared to the plates but has similar 
density to the plates, it is safe to say that the shield will 
only have a thin-shell when the plates have thin-shells. 
Assuming the plates do have thin-shells, we define m s to 
be the mass the chameleon would have deep inside the 
shield if the shield has a thin-shell i.e. m s — m^tfrs) 
where V'(<p s ) = —Pps/Mp\- Since the shield is sand- 
wiched between the two plates, the thin-shell condition 
for the shield is simply m s d s > 1. When the shield has a 
thin-shell, its presence attenuates the chameleonic force 
and torque on the detector due to the attractor by a fac- 
tor of exp(— m s d s ). Since exp(— m s d s ) w 1 in the absence 
of a thin-shelled shield, we can take account of the shield, 
thin-shell or not, by changing the definition of thus: 



Tj,(d) w a T e~ 



EM 

A 



(41) 



This expression provides a very good approximation for 
theories in which the precise value of rh is unimportant 
(e.g. ones with rribrh 1) and an order of magnitude 
estimate otherwise. 



C. Inverse Square Law Constraints 

The 2006 Eot-Wash experiments requires that 
\T (j> {d= 55 /im) | < 0.87 x 10" 17 Nm, 
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with 95% confidence. We define T^(g? = 55/im) = o,tA t 
and find that the above bounds correspond to: 



A T < 0.89 x lCT 12 GeV, 



(42) 



Importantly this is smaller than the energy scale as- 
sociated with dark energy: A c i c = 2.4 x 10~ 12 GeV; 

Pde = A^ . 

Using our expression Eq. ([4Tj) for the chameleonic 
torque, we find that the constraints we must apply are 
as follows: 



55 fim 



A 



ds < 7.0 x 10" d ' GeV d .(43) 



Assuming that m b <C m c where m c is the chameleon mass 
inside the plates, it follows that F^/A has the following 
form: 



A 
where 



In 



Ro 



R b 
Rq 



1 



1 



m 2 £ Q 



So 



2C{m b r p )m b 



■m c 



mlC{m b r v ) 



AC(m b r p )m b 



(47) 



When C(m b r p )m b /m c « 1 we have So ~ 2C 2 m 2 /m 2 
and so the last term in Eq. (147)) is: 



V. APPLICATION OF CONSTRAINTS TO f(R) 
THEORIES 

A. Chameleonic force 

The chameleonic force per unit area between two par- 
allel plates is given by Eq. (|34p. To prevent large de- 
viations from general relativity occurring over solar sys- 
tem, and smaller, scales, one must require that f(R) ~ 
R + h(R) where \h'(R)\ < 1 and \h(R)/R\ < 1. In this 
case the expression for F^/A becomes: 

^ « ^ i(Ro - Rb) h'(Ro) + (h(R b ) - h(Ro)) + To] ■ 

where for R c 3> R b 

R 2 c h"{R b ) 
4C(m b r p ) 

£ = 1 + 2C(m b r p )V (R c , R b ) - + ACV , (45) 
h(Rc) - h{R b ) - {R c - R b )h'{R b ) 



V 



h"{R b )R 2 c 



(46) 



where r p is the radius of the parallel plate(s). We shall 
now consider several potential forms for h(R). 



B. Logarithmic potentials 

We begin by considering a simple chameleon gravity 
model that was recently suggested in Ref. [37f for a gen- 
eral /?. The theory, when written as a chameleon theory, 
would have a potential V(<p) = Vo — A 4 , ln(0/Mpi), it 
was suggested that this would result in an experimen- 
tally viable and cosmologically interesting dark energy 
model, where Ag/M^ ~ 0(Hg) [13 for j3 < 1/4 V3. We 
will analyse the same model in the f(R) setting where 
(3 = 1/ %/6 > 1 /4V3 and show that local tests already 
lead to difficulties, see also [38l |. 

On laboratory scales we would have f(R) ~ 1 + h(R) 
and so we find: 



h(R) 



2A| 



A n ,.,-2 



11)- 



-S a w A 4 C(m b r p ). 



2mgC(m 6 r p ) 

Alternatively if C(m b R)m b /m c > 1 we would have: 



2m 2 C(m b r p ) 



So 



A 4 - 



If C(m b r p ) -C 1 i.e. m b r p 3> 1, then it is clear that this 
last term is always small compared to the other terms, 
however if m b r p -C 1, then the last term will dominate 
the expression for the force. 

The chameleon mass for a given R in this set-up is 



m^R) 



In between the two plates, 
d 2 ^ 



MpiR 

: 71a 2 ' 

satisfies 311: 



dz 2 



and (f>o is defined to be value of <p midway between the 
two plates (i.e. a distance d/2 from either plate), where 
by symmetry d(f>/dz = 0. Integrating the above equation 
we therefore have: 



^) =2(v(q>)~v(M-v4<j> b )(t~M). 



Integrating this again and defining <f> s ~ 0((j> c ) to be the 
value of <p on the surface of the plates, we have: 



dx 



V2 4 v^FW^MFW' 

Following Ref. [3l| , when m b <§; mo <C m c we have that 
4> s ~ 0(<p c ) < (j> < 4>b and so V '{$)- V{<j) )-V^{<j) b ) {<j>- 
O ) a V{<j>) - V{4>o) = A 4 o M0o/& and so: 



d 



A 2 



da; 



o Jo y/ln(l/x) 



Noting that mo = Aq/0o and evaluating the integral we 
find: 



mod = 



M Pi R d 
VGA 2 



2lt . 
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Now mj « mo « m c is clearly equivalent to R c ^> Rq^> 
Rb, and in these cases we therefore have: 



Ro = 



Mpid 



(48) 



It follows that, irrespective of the value of £q, F^/A 
drops off more slowly than 1/d for all m^d <C 1. In 
the latest version of the Eot-Wash experiment [Hj], the 
plate radius, r p , is 3.5 cm, and the smallest hole radius 
is 1.6 mm. The pressure of the laboratory vacuum is 
10~ 6 torr which corresponds to a background density of 
p b = 1.6 x 10~ 9 kgnr 3 = 6.7 x 10- 30 GeV 4 . 

Now if the vacuum region is large enough then R b = 
Pb/Mp! and so nib = fhb = Pb/vwVnMpi. However, it 
was shown in Refs. [28l . l30l . [31| that if the vacuum re- 
gion only have length scale L vac and fhb <C l/L vac , then 
generically mj, ~ C(l/L vac ). For the moment we only 
assume that r p /L vac <C 1. 

We therefore find that for d = 55 pm we have rribd < 1 
for all A > 5.6 x 10~ 19 GeV, m b r h < 1 for A > 3.0 x 
10~ 18 GeV and m b r p < 1 for A > 1.4 x 10~ 17 GeV. 
The suppression factor due to the electrostatic shield is 
exp(— m s d s ) where: 

, gPshiekA non { 10- 12 GeV \ 2 

where we have used d s = 10pm and p s — 8.3 gem . 
Thus, whenever m s d s < 1, we are therefore firmly in the 
uibT-p iribr^ <§; 1 region and hence C(mf,r p ) 3> 1. Thus 
C(m b r p ) ps l/2m 2 r 2 and: 



A 



— In(mfcti/v27r) + mbd/v2it — 1 



From Eq. (|4ip , in the absence of the electrostatic shield, 
the chameleonic torque for 3> d would be: 



T<j, ss drAgr/i 



In 



2tt 



(49) 



Here m c is the chameleon mass inside the plates which 
have density p c pa 10.2gcm~ 3 . 

We note that the requirement that the plates have 
a thin-shell constrains the value of mj, and it is im- 
portant to check that this constraint holds. Conser- 
vatively, the thin-shell constraints for the plate require 
{3(<frb — <f>c)/Mp\ < $at/3 where $tv is the Newtonian 
gravitational potential of the whole experiment at the 
surface of the plate. Since the geometry of the experi- 
ment is complicated, we do not calculate $jy Instead, we 
estimate <S> N /3 < 10~ 26 , and so fa - C < 7 x 10~ 8 GeV. 
Given that p c 3> Pb, we take (f> c -C <f>b an d then from 
mb = Aq /</>{, we must have 



1/mj < 14i 



10~ 12 GeV 

An 



(50) 



If this condition does not hold, then the plates would 
not have thin-shells and the Eot-Wash data would au- 
tomatically rule out the theory. The experiment takes 
place inside a vacuum chamber with smallest dimension 
ivac = 0.2 m [33| . We assume that the walls of the vac- 
uum chamber have thin-shells. Approximating the walls 
of the vacuum chamber perpendicular to the shortest di- 
mension as being parallel plates, we use Eq. (|4"8|) above 
to tell us that when the background density of matter in 
the vacuum chamber is very small, we have in the centre 
of the chamber: 



Rh Rvac 



Mp\L vac 

This formula holds as long as Rb = pb/Mpj < -R vac - In 
the opposite limit we just have Rb = R b - In all cases we 
have Rb > R v &c and so: 



m b > 



/2vr 



Mp\R vac _ 



12m' 



(51) 



and so condition (jST))) is always satisfied for Ao < 1.2 x 
10~ n GeV. 

Given Eq. (f5Tj) . we find that for the allowed values of 
Ao we always have: 



C(m b r p )m b 



1 



2m b m c r 2 



« 1, 



where we have used m c = p c / \/l!>Mp\K^. Thus: 

£ ps 2C 2 m 2 /m 2 = (77i c m b r p )~ 4 /2. 

It follows that, in the absence of the electrostatic shield, 
the chameleonic torque for r% d is 



a.TA r h 
1 

+ 



In (y2Trm b rh^ 



m b r h 



2tt 



(52) 



2m 2 r 2 



Including the suppression factor due to the electrostatic 
shield, which is exp(— m s d s ), we therefore find the fol- 
lowing constraint on yo = Aq/(10~ 12 GeV): 



which gives yo < 0.37 and so: 

A < 3.7 x 10~ 13 GeV. 



(53) 



Cosmologically, the mass of the scalar field at the min- 
imum of its potential is given by 



3 £l m M p iH 
2 ~AT~ 



and the value of (j) at this minimum is given by: 



xmin 
cos 



A 4 



Mi 



pi 



3Vl m H 2 M 2 { 
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Now H = 2.1hx 10" 42 GeV and from WMAP 39]: Q r , 
0.127/i~ 2 and h = 0.73. We find that 

m cos /H = 

Va 



min 

cos 



pi 



0-06y 4 



(54) 
(55) 



and so the Eot-Wash constraint on Ao gives: 
m cos /H > 8. 

At its minimum then, the <f>- field is still heavy today. This 
should be contrasted with the requirement obtained in 
[37l | that the mass of the 0-field should be small compared 
to the Hubble rate in order to drive acceleration. Here we 
find that local tests and the thin shell requirement impose 
that the mass of the 0-field at the cosmological minimum 
is so large that the field must sit there on cosmological 
scales. It is easily checked that m cos /H is a decreasing 
function of time, and so in the past (f> was heavier still 
relative to H. Therefore 4> will have remained stuck close 
to the minimum of its effective evolution throughout the 
matter era. Additionally the Eot-Wash constraint on A 
implies that: 



< 0.001 



pi 



We have considered the potential V(<p) = Vb — 
Aq1ii(<^/Mpi). The constraint on Ao implies that 
A^/SM^H 2 < 0.00026 and so if V{<f>) is to be the source 
of dark energy and there is to be a realistic amount of it 
today, we would need V > O(1000)A 4 . Note that this is 
very different from the original scenario envisaged in ref. 
[3^ . where Vq ~ A 4 , so that the whole potential could 
be written in the form V = — Aoln(0/M) where M ~ 
0(Mp\). The Eot-Wash constraint on Ao therefore rules 
out a scenario where Vq ~ Aq for (3 = 1/V&, confirming 
the cosmological obstruction noted when (5 > l/4\/3. If 
we moved away from f(R) theories and allowed for dif- 
ferent couplings, we would find similar constraints on A 
for other 0(1) values of the coupling f3. 

Relaxing the constraint Vq ~ A 4 and allowing much 
smaller values of Ao, it should also be noted that the 
conservative thin-shell constraint for a test mass with 
density » Aq on the cosmological value of <j> (as derived 
in Section IIIip actually provides a stronger constraint on 
the cosmological value of the field today and as such gives 
a tighter bound on A . Specifically Eq. (f31~j) implies: 



Mpi 



< 10~ 4 & A < 1.8 x 10^ GeV. (56) 



This leads to the following constraint on the mass of 
at its minimum cosmologically: 



Since m cos /H ^> 1, <fi lies close to its cosmological mini- 
mum and so, in the Jordan frame, by Eq. (|18[) : 

-2$ 3 U$ = ^P* 2 ^ - -Pmatter, 

where $ = e -W/ M Pi and so $ « 1. To leading order 
with p = lndwe have <p p « 3</>. Therefore 



$ Mpi ~ Mpi 

and to the same order 



6Ag 

/^matter 



4> 



pp 



18/30 



pi 



I8A0 _ 

/^matter 



S /H > 35. 



(57) 



We define 9 = A 4 /p mat ter and then using Eq. (|2"T|) we 
arrive at 

(1 + w cS )il dc « 3fl-i ' ' 

l-36> 

where 6*o is the value of 8 at the present time. Assuming 
that the Universe is flat (k = 0) and taking — 0.76, 
we find today when t = t^: 

l + w cS {t Q ) « -4.3200- 

Notice that the effective equation of state is below -1, 
this is a consequence of the scalar-tensor character of the 
chameleon model. 

The Eot-Wash constraint on Ao gives 

|l + w ff| < 0.0085, 

while the thin-shell constraint on Ao gives 

|l + w cff | < 10- 4 , 

which is in line with our expectations from Section ED 
Whilst the thin-shell constraint on the cosmology is much 
stronger than the Eot-Wash bound, the cosmological con- 
straint makes a number assumptions above the nature 
of inhomogeneities in the Universe, in particular about 
their scale at the present time. One could presumably 
argue that the cosmological constraint could be relaxed. 
The same line of argument cannot be used for the Eot- 
Wash constraint, and as such represents a strong con- 
straint on the magnitude of deviations from ACDM in 
this model. As such, the model cannot be distinguished 
from a ACDM model at the background level. At the 
perturbative level, the situation is very different as the 
bound (|57[) implies that density contrast would have an 
anomalous growth on scales lower than 100/i _1 Mpc. This 
may be testable in the near future with next generation 
red-shift surveys [25|. 

The version of the logarithmic potential f(R) theory 
suggested in Ref. [37| required A 4 . > p m attcr today i.e. 
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Ao > 1.7 x 10~ 12 GeV. Even from a conservative point 
of view (and indeed for any (3 where 2f3 2 ~ 0(1)), such 
a value of A would produce a torque in the Eot-Wash 
experiment that is almost 100 times larger than the 95%- 
confidence level upper bound. The scenario suggested in 
Ref. [37j is therefore strongly ruled out by local tests of 
gravity. 



C. Power-law form 

In many cases [H|, IH, one finds that for R > H 2 , 
where Hq is the Hubble constant today, h(R) has a power 
law form i.e.: 



h(R) 



P 



+ 1 \R 



P+i 



(59) 



for some p ^ and some constant R > 0. For a 
chameleon mechanism to exist we need V < 0, V" > 
and V" < 0, and so must require p < 1. Relative di- 
vergences from GR such as those parametrized in the 
PPN formalism or measured by observing the motions 
of planets would scale as h(R)/R and h'(R) or by the 
ratio of any variable component of the effective cosmo- 
logical constant to the local matter density . However, 
the Eot-Wash test probes changes in V(cf>), which scales 
as h(R) and Rh'(R), although they are only sensitive to 
this when the chameleon mass in the background, which 
scales as l/h"(R), is not too large. 

In theories with < p < 1, both h(R)/R and V(4>) 
would be largest for large values of R. These theories 
would therefore diverge most markedly from General Rel- 
ativity in the UV (i.e. large R) regime. Increasing R 
would make both h(R)/R and h"(R) smaller, and so ulti- 
mately one could ensure compatibility with all laboratory 
tests by making R very large. Provided h" is not small, 
however, the changes in V(<f) that could be detected by 
the Eot-Wash experiment would increase. 

If -1 < p < then h(R)/R and h'(R) are largest 
in the IR regime where R is small. However V(<f>) still 
increases with R, and since R increases as the separa- 
tion of the plates in the Eot-Wash experiment is de- 
creased, the smaller the separations the stronger the po- 
tentially detectable signal would be. Ultimately compat- 
ibility with all local tests could be ensured by making R 
small enough. Additionally in all theories where p > — 1, 
F$/A would be dominated by the ^-dependent (i.e. Rq 
dependent) terms and only weakly depend on Rb when 
m b d <C 1. 

Finally, in theories withp < —1 both V(<f) and h(R) /R 
would decrease with R. This would mean that F^/A 
would only very weakly depend on d and generally be 
much smaller in a given set-up than for the other classes 
of theories. Again compatibility with all local tests could 
be ensured by making B small enough. 

The —1 < p < theories are the most testable type 
of theory as they would result in deviations from GR in 



both the UV and IR regimes. In the UV regime there 
would be potentially detectable fifth-forces between par- 
allel plates, and in the IR regime the ratio of the density 
dependent part of the effective cosmological constant to 
the ambient matter density would increase cosmologically 
at late times as the ambient density decreased. 
In all of these theories: 



A 



M^Rop 2 (B, 



2(p + l) V R 
M^Rbp (RbV _ M^RbP fRo' '• 
2{p + l)\R) 2 \R 

M 2 x p 2 R 2 c (R b \\ 



(60) 



8C{m b r p )R b \ R 
where £q is given in terms of C and T>q by Eq. (I45|) and 



Vn = 



1 



p( P +l) 



1-p 



Rb 

R, 



Note that the last term in Eq. ()60|) is independent of the 
plate separation d and vanishes in the limit m b r v — > oo. 



1. Relationship to chameleon theories 

Converting these theories to chameleon theories we 
have for h' <SC 1: 

- 2 wr h{R)=p {R) ' 

and so R cx (p 1 ^ and h(R), Rh'(R) cx ^p +1 ^p. It follows 
that 

VW = const + (Z^fy. 
W 2(p + 1) V M PlP J 

and so defining n = — (p + l)/p, we see that, neglecting 
the constant term in the potential, |V(</>)| cx |0| _n - In 
the context of chameleon theories these potentials have 
been studied in great detail 0, H^, [2t| [3(1 HH , and so 
we are able to apply a raft of results to the analysis of 
these theories. 

In appendix [Bj we show the mass of the chameleon 
field for m c 3> tuq 3> mb (w c is the chameleon mass 
deep inside the plates and nib is the chameleon mass in 
the background) is given by mod = a p where: 



p 2 B 



1+p 



P 



-B 



2' (1+p) 
1 (1-p) 



P<-1, 
1 <p< 1. 



Using mo 
Bq <C R c '- 



1+p \2 2(l+p), 
l/3/i"(i?o) we therefore have when Rb <C 



2„2 \ 1-p 



3p 2 a 
Rd 2 



(61) 
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For what follows we also define: 



Constraints on Power-Law f(R) theories: f(R) = R+h(R) 



K p = (3p 2 a 2 p ) lt » . 



(62) 

Using rn ~ l/3/i"(-Ro) we therefore have when <C 
Ro <C Rc'- 
XJsmg the relationship between Rq and d derived above 
when R c > R > R b Eq. (|6"0|) becomes: 



1 \ "-»■ 

RePJ 



M^KppiR 
2(p+l) 

8C(TO h r p )i? b \ R J 



G a 



(63) 



£ - 



where 



G 



( mt,d 
p \a p 



= 1 + - 



1 / rribd 



2(l + p) 
1-P 



(64) 



p + 1 ( nibd\ 1_p 
P 



We note that G p = when m^d/ a p — 1 which corre- 
sponds to Ro — Rb- We now consider the integral: 



/(d,rh) 



.4 



-ds 



The approximation used to calculate Rq (d) breaks down 
when R « i? b , which corresponds to rribdja v ~ 1. In 
the case mbrh/a p > 1 we cannot simply use Eq. (|6"3"|) to 
calculate /(ci, rh) as we must integrate over values of d for 
which Eq. (|63|) is not valid. However we should be able to 
trust Eq. (|63|) for smaller values of d. For m,bd/a p 3> 1 
we expect an exponential drop-off in the force, just as 
one would find in a Yukawa theory at distances larger 
than the inverse mass of the scalar field. We therefore 
do not expect the dominant contribution to I(d, r^) to 
come from values of d < a p /mb- We also note that if 
rribrh/ap ^ 1 then TOi,r p > 1 as ^ < r p and as such the 



second term in Eq. (|63[) is negligible. The first term in 
Eq. (|63p vanishes when nibd/a p = 1, and since we do not 
expect a significant contribution to the integral to come 
from larger separations, we evaluate I(d, r^) by using Eq. 
(|6"3"|) for F < p(s)/A but if TO(,rh/a p we cut the integral off 
at a separation a p /mh. 

Thus we define x(d) = m,bd/a p and x max = 
imn(mbrija p , 1) and find: 



I(d,r h ) 



p 2 M 2 1 K p R 1 / 2 



TOfc 



RV*Or 



l + 3p 
1-P 



(65) 



x {[H p (x(d)) - H p (x mlix _)) 



4R 2 b C(m b r p ) 



■So 



where 



H p (x) 



1 



1 



(1+p) 
1 — p 



1 



- p i+3p a; 

— x i-p 

3p p 



l-3p 

a; !-p 




Allowed Region 


Eot-Wash bound 

(thin shells assumed) 
- - - Cosmological thin shell bound 
m D » 1 




c p 



Slope of h(R): p 
Constraints on Power-Law f(R) theories: f(R) = R+h(R) 




Allowed Region 



Eot-Wash bound 

(thin shells assumed) 
■ - Cosmological thin shell bound 
, i m D » 1 



-0.E 



-0.6 



-0.4 



-0.2 



Slope ofh(R):p 
Constraints on Power-Law f(R) theories: f(R) = R+h(R) 



Eot-Wash bound 
(thin shells assumed) 
- Cosmological thin shell bound 




Excluded Region 



Slope of h(R): p 

FIG. 1: Eot-Wash constraints (thick solid blue line) on 
f (R) gravity theories with f(R) = R + h(R) where h(R) = 
R(R/R) P+1 ; R = A(5/Mi, and -5 < p < -1, -1 < p < 
and < p < 1. For this constraint we have assumed that the 
test bodies have thin-shells (which is necessary to avoid local 
tests). We have also shown : (1) the cosmological thin shell 
constraint (thick red dashed line) for test bodies in the labo- 
ratory derived in Section Mil (2) the naive constraint (thick 
black dotted line) one could derive by simply requiring that, 
inside the test bodies, the mass of the chameleon at the min- 
imum of its effective potential, m c is large compared with the 
length scale of the body, D p . This was the constraint con- 
sidered in Ref. [34} . For all such theories we see that the 
correctly evaluated constraint provided by the Eot-Wash ex- 
periment [32] is stronger than both this nai've constraint and 
the cosmological thin-shell bound for all for p > —1. The 
m c Dp ^> 1 constraint never provides the strongest constraint. 
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With this formula we are able to evaluate the Eot-Wash 
constraint for all theories with h(R) oc R p+1 . We do this 
further below. However, we discuss first the cosmological 
thin-shell constraint on these theories. 



2. Cosmological Constraints 

On cosmological scales, the field <fi is stuck at the min- 
imum of the effective potential provided rrA/H 2 ^> 1 
which becomes: 



4» dc + n m 

Rh"(R) 

If h(R) oc R p+1 this becomes: 

4^ dc + » m 
\ph'(R)\ 



» 1. 



> 1. 



The cosmological thin-shell constraint requires that: 
I3\A4>\ 



Mpi 



< 10" 



where Acf> is the difference between the value of 4> cosmo- 
logical and the value of 4> at the minimum of the effective 
potential in a region with density 0(1) g cm~ 3 . This gen- 
erally implies that cosmologically f3\cf>\/Mp\ < 1CP 4 and 

k\ h '( R )\ < 10 ~ 4 - lt is clear then that E( l- © holds 
provided p x 10~ 4 <C 1 and so for 0(1) values of p, we 
are always in the region where m 2 ^/ H 2 3> 1 and cf) lies 
close to the minimum of its effective potential. 

To leading order we take ~ — /3p m atter/-Mpi and, 
defining p = In a, where a is the FRW scale factor in the 
Jordan frame and $ = e ~W/ M Pi we nn d that: 



3. Collected Constraints 

We will consider now how the Eot-Wash data, when 
thin-shells are assumed, constrains the properties of 
power-law f(R) theories. It should be stressed that in 
the absence of thin-shells, the Eot-Wash would automat- 
ically rule out these theories. 

Defining R = A 4 /Afpj we have plotted the Eot-Wash 
constraints on Ao for — 5 < p < 1 in Fig. [1] as a thick 
(blue) solid line. The cosmological thin-shell constraint 
is shown as a thick (red) dashed line. For theories with 
< p < 1 we find a lower bound on Ao and for theories 
with p < we recover an upper bound. We also show, 
as a thick (black) dotted line, the naive constraint on 
the parameters that one would find by simply requiring 
that the chameleon mass at the minimum of the effective 
inside the plate, m c , is large compared to the plate thick- 
(66) ness D p = 0.997 mm. It is a commonplace assumption in 
the literature (see e.g. @,[33|) that assuming m c D p ^> 1 
(where more generally D p would be the length scale of 
the test body) is enough to satisfy local tests of grav- 
ity. It is clear from the plots that this naive bound never 
provides the tightest constraint on the parameters of the 
theory 

The constraints on R constrain the equation of state 
parameter of the dark energy described by the f(R) the- 
ory. Taking fi m = 0.23 today, we plot the collected con- 
straints on the effective Jordan frame equation of state 
parameter (as defined in Section [TTT]) in Fig. [21 We see 
that at the current epoch |1 + Wdc| < 10~ 4 for all 0(1) 
values of p with the largest values occurring for p < — 1, 
and hence the late time cosmology produced by any vi- 
able theory would be observationally indistinguishable 
from that described by the standard ACDM model. 



3ri m H 2 



where /o = £l m / (^ 
<I> 



f « 9f Rh" 



-3f Rh"(R) < 1, 



f Rti"(R) 



h"{R) 



Today from Eq. ([27]), with fi^f = fl m , n 
1 — fl m we have /o « £l m / (4 — 3f2 m ): 



off 
de 



(67) 

(68) 
f2 de « 



(1 + wf c )n dc « 3f Rh"{R) 



4 foRh'"(R) n„ 
3 h"(R) 2 



and so for f2 m = 0.24 and h(R) oc R p+1 we have: 
|1 + wf e \ftf e w 0.32\ph'(R)\ |1 + 0.050(p- 1) 



(69) 



The cosmological thin-shell constraint ensures that cos- 
mologically \h'(R)\ < 10~ 4 today and so: 

\l + wf e |Of <3.2|p||l + 0.050(p-l)| x 1(T 5 . 



VI. CONCLUSIONS 

In recent years, modifications of General Relativity 
have been suggested as a possible explanation for the ob- 
served accelerated expansion of the universe. A popular 
class of models are the so-called f(R) theories. While 
cosmologically viable theories can be found, local con- 
straints on such theories have to be worked out, since 
the gravitational sector is modified, which could result in 
unacceptable deviations from Newton's law of gravity. 

In this paper we have constrained f(R) theories, us- 
ing the well known equivalence between these and scalar- 
tensor theories. For an f(R) theory to be consistent with 
both cosmology and local gravity experiments, the equiv- 
alent scalar-tensor theory must be a chameleon field the- 
ory. We have shown that the requirement of the thin- 
shell mechanism at work in Eot-Wash experiments re- 
sults in an equation of state for dark energy very near to 
that of a cosmological constant. Thus, viable f(R) mod- 
els (those which are compatible with local experiments) 
behave on cosmological scales similarly to the standard 
ACDM model and deviations are expected only on very 
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Combined constraints on Dark Energy Equation of State 



- Upper bound on | 1 + w 



_7 _e _5 _4 _3 _2 -1 
Slope ofh(R):p 



Combined Constraints on Dark Energy Equation of State 




Slope of h(R):p 
Combined Constraints on Dark Energy Equation of State 



« 10- 
E 



- Upper bound on | 1 + w 



Slope of h(R): p 

FIG. 2: Combined Eot-Wash constraints on the effective dark 
energy equation of state parameter produced by f(R) gravity 
theories with f(R) = R + h(R) where h[R) = R(R/R) P+1 . 
R = Ao/Mp! and -5 < p < -1, -1 < p < and < p < 1. 
These constraints have been derived by requiring both that 
the Eot-Wash test masses have thin-shells and by requiring 
that the chameleonic torque produced between the two thin- 
shelled test masses is small enough to have avoided detection 
to date. We see that in all cases we have |1 + < 10~ 4 
today. As a result, the late time cosmology of any viable the- 
ory would be virtually indistinguishable from that described 
by the ACDM model 



[37l | for a fixed coupling (3=1/ \/6) as well as power-law 
potentials (such as those presented in [H, HH). The for- 
mer are ruled out by local gravitational tests, while there 
is still room for the latter models. 

To conclude, while on cosmological scales viable f(R) 
theories behave like ACDM, deviations are expected on 
scales which could be large enough to be within the reach 
of next generation galaxy surveys [2fj| . Hopefully, future 
measurements of the dark matter distribution on those 
scales can be used to find such deviations from the stan- 
dard ACDM model. For this, a detailed understanding of 
galaxy formation is necessary, including an understand- 
ing of both the dynamics of baryons as well as that of 
dark matter in ACDM and /(i?) /chameleon theories. 
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APPENDIX A: THE FORCE BETWEEN TWO 
PLATES 

In previous works (28l . l3l| , the chameleonic force per 
unit area between two parallel plates was calculated and 
found to be: 

^ = y'(to)-y'(&) + |gfa ) -&), 

where 4>$ depends on d. We assume that both plates have 
radius R and thickness D and that D <C R- However as 
these calculations treat the plates as being infinite, they 
required for consistency that: 

• either, for an isolated plate, <f> — > <pb at a distance, 
d <C R, from the plate so that the infinite plate 
approximation was still valid, 

• or, the precise value of <pb was not important when 
the plates were separated by a distance d <C R. 
This means that provided V'(<f>o)/V ((/)*) <C 1, one 
could replace <pb in the above expression by <jf with- 
out altering the prediction for F^/A greatly. In 
these cases, the behaviour of <f> far from the plates, 
where the infinite plate approximation is invalid, 
would be unimportant. 



small (sub-galactic) scales. The expected deviations from 
the cosmological constant equation of state w = — 1 now 
in viable f(R) theories are unmeasurably small (at least 
with current technologies). As examples, we have stud- 
ied f(R) theories with logarithmic potentials (based on 



These approximations held for all of the chameleon the- 
ories considered in Rcfs. [28], [3l|, however in this work 
we consider a wider range of theories, and it is often the 
case that both of these assumptions fail to hold. In this 
appendix, we therefore derive an improved version of the 
force formula. 
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Outside of a body in a region where the background 
density is p b , the chameleon field obeys: 



V 2 = V^cj)) + 



Ah 



pi 



Consider this equation near one of the circular surfaces 
of a cylindrical plate, of uniform density, and with radius 
R and thickness D <C R. Defining z to be the distance 
from one of the circular surfaces of the plate we have: 



d 2 4> 
dz 2 



ld_ 

r dr V dr 



(Al) 



We begin by consider the case where only one plate is 
present. Here 4> — ► <fib, dip/dz — > as z — > oo. Integrating 
Eq. (|Aip with these boundary conditions give: 



+ / i£(r£)£d,. 

r or \ or I dz 



(A2) 



We solve this approximately by assuming that for z < 
z* the z-dependence of the r-derivative terms is weak 
compared to that of the potential terms, and that for 
z > z* , the non-linear terms in the potential, i.e. terms 
that depend on 3rd or higher derivatives of V , are sub- 
dominant. In z > z* , we have 6^6 where: 



~\ 2 




(A3) 



ld_ I d± \ d± dz 
r dr \ dr I dz 



or equivalcntly: 

dz 2 



ld_ 

r dr 



dr 



Assuming that the plate is thin (D <C R), and solving 
Eq. (|A3jl . we find that for the z > z* and along r = 0: 

- <t> b oc (e- m " z - e- m6r ) , 

and here z is the distance from the plate surface. It 
follows that for z <C R we have: 



i d I d<f>\ dc^ 2 

r or \ dr dz 



,(e m » R - 1/2) 



(A4) 

In z < z* , we have assumed that the z-dependence of the 
r-gradient terms is relativity weak. We therefore approx- 
imate the r gradient terms in Eq. (|A2|1 using the z > z* 
solution i.e. we approximate them using Eq. (|A4|) with 
<f> — ► 0. For z i? we then have 

2 



2 I fl: 



-ml((j)- (f> b f 



(e mbH - 1/2) 

( e m,b-R _ 1)2 



(A5) 



The above equation also holds approximately in the 
z > z* region, provided z <C i?, and so provides an 
approximation to the evolution of (f> everywhere when 
z -C i?. In particular we see that on the surface of the 
plate, at z = 0, where say = S : 



y(0 s )-y(</> 6 ) + v (0 b )(0 s -0 b ) 

(e mi>i? - 1/2) 



0' 



( e m 6 ii _ 1)S 



(A6) 



We assume that the plate has a thin-shell, so that deep 
inside it </> — > C where: 



f3 P c 

Mpi ' 



Provided the shell is thin, we can treat the system as 
being essentially 1 dimensional [28|, [31| and so: 



V(0)-V(^c)-^(0o)W-0c). (A7) 



Thus by evaluating and equating the left hand sides of 
Eqs. (|A6[) and (|A7[) at the surface we find: 



-v^{fhW.-<k) 

+m 2 b C(m b R)$ s - fo) 2 = 0, 



where 



Thus: 



C(m b R) 



(e 



m b R 



1/2) 



^ e m b R _ 1)2 
(^(0fc)~ Wc)) 

2m 2 C{m b R) 



where 

D(<p c ,c/) b ,m b R) = 



-y/l + 4C(m b R)D((/> c ,(l> b ) (A8) 



(^(<M~^W>c)) 2 



We now consider the force between two parallel plates. 
This derivations make uses of results found in Refs. 28] 
and [31] , and proceeds along roughly similar lines. 

In between two parallel plates with radius R and with 
separation d <SC R in the z-direction, the chameleon field 
obeys: 



d 2 cf> 
dz 2 



(A9) 



For simplicity we treat the plates as having the same com- 
position. This assumption was dropped in Ref . [281 ] , how- 
ever, it was also shown there that for most purposes the 
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assumption provides an excellent approximation. This is 
because the chameleonic force generally exhibits very lit- 
tle composition dependence [23]. We define z = to be 
the surface of one of the plates, and z = d to be the facing 
surface of the second plate. The system is symmetric and 
so &4>/&z = at z = d/2. We define <fi(z = d/2) = 4>o(d). 
A formulae for 4>o(d) have been provided in Ref. [3l| . 
Integrating Eq. (|A9|) we have: 



v(<M-va<MU>-<M- (Aio) 



Following Ref. [31|, inside either plate, Eq. (|A7|1 holds. 
By equating both Eqs. (|A10j) and (|A7|) at the surface of 
one of the plates, where <f> = <f> s say, we find: 

, _ [v{<t> c ) - v{<h) + vA<M<fo - vA<f>c)<Pc} 
vA<k)-vA<i>c) 

In Ref. [3l| it was shown that the attractive 
chameleonic force unit area between two thin-shelled 
plates is given by — V^{(j) c ){4> s — (f> s ), and if, as is usually 
the case, the plates are much denser than their environ- 
ment so that V^(4> c ) /V^{4>b) = Pel Pb 3> 1, we have: 

^ - V{4>o)-V^ b )~Vl{4> Q -4> b )+mlC{m b R)(4> s -ct> b f. 



(All) 

ami 



This coincides with the formulae found in Refs. [26 
when m b R 1 => C(m b R) ps 0, or more generally when- 
ever the final term is small compared with the other 
terms, which is when C(m b R)D((f> Cl (f> b ) <C 1. 4> s — <t> b 
is given by Eq. (jA8|) . 

When m b R < 1, we have C{m b R) « \/2m 2 h R 2 > 1. 
If this is the case we also have C(m b R)D(cf> c , tp b ) ^> 1: 

m 2 C{m b R)$ s -fo) 2 « F(^)-F(0 c )-F (^ c )(^-^ c ), 

and then in this limit Eq. (|A11[) becomes: 

Fa, 



A 



= V(M - V{<t> c ) - V&fa - <f> c ) - v b \<Po - fa). 



Thus integrating the 4> equation we find: 



In this work we have considered power law potentials 

p+i 

where V cx ep(—e<p/(p + 1)) ? and where e = sgn(p(p + 
1)) and p < 1. For these potentials we have m| = = 
(p + I) / p 2 V ((/))/ 4> 2 . Thus, defining Y — (j>/4>o and m = 
in4>(4>o) we have: 



1 (dY 

2 \~dz~ 



2 2 

p TUq 

(p+1) 



p+1 



(B2) 



-1 



p+1 



(y-i 



Now 



R oc cf> 



l/p 



and so defining i?, = R((j>o) we have ((/>b/</>o) 1 ' /p = Rb/Ro- 
Thus when R b /Ro <C 1, the last term in Eq. (|B2j) is very 
small and can be dropped. Working in this limit we have: 



dY 
~dz~ 



2 2 2 



y— - i 



Now on the surface of the plates 4> ~ 0(<p c ) where <^) c is 
the value of 4> inside the body, and assuming R c ^> Rq 
i.e. m c 3> mo, we can we treat (4>/4'o) p as becoming very 
large as z — ► (i.e. as we approach the surface of the 
plate). We define X = Y~ p = (<j>/(j)oy p and then in the 
limit m b <C mo <C m c we have: 



dA 

dz 



2m 2 X 2+2 /P r _* 
X > 



(p+1) 



APPENDIX B: CHAMELEON MASS BETWEEN 
TWO PLATES 

In this appendix we generalize the calculation of the 
chameleon mass between two parallel plates, as per- 
formed in Refs. [H, [31]] > to include the wider range of 
chameleon theories considered here. 

In between two parallel plates with, say, a circular cross 
section, in the x — y plane, with radius r p , and a sepa- 
ration, d, in the z-direction where d <C r p , Eq. ((6]) for <\> 

simplifies to be essentially one dimensional □ 



dz 2 ■ 



d 2 

dz 



(Bl) 



where <fi b is the background value of <; 
to be the value of <fi when d<fi/dz = 0, 
midway between the two plates when z 



i>. We define <f>o 
which will occur 
= d/2. 



Integrating this equation and using X = at z 
X = 1 at z = d/2 we have for p < — 1: 



mod 



P 



;Bll 



P) 



,/2\p~tt\ \p+i\ V2'(i+p);' 

which simplifies to: 



mod 



, |1+P 
If — 1 < p < 1 then we find 



2 Ml- p 



2'(1+P) 



m O d =\lJT— — \ B \\i- 



l+p| V 2 ' 2 ( 1 +P) 



(B3) 



(B4) 
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